DISCRETE SCALE INVARIANCE IN TURBULENCE? 



00 
PL, 



43 

o 



CO 

c3 



D. SORNETTE 1 * 2 
1 LPMC, CNRS UMR6622 and Universite de Nice-Sophia Antipolis 
B.P. 71, 06108 NICE Cedex 2, France 
2 IGPP and Department of ESS, UCLA 
Los Angeles, CA 90095-1567 
(February 1, 2008) 

Based on theoretical argument and experimental evidence, we conjecture that structure functions of 
turbulent times series exhibit log-periodic modulations decorating their power law dependence. In 
order to provide ironclad experimental evidence, we stress the need for novel methods of averaging 
and propose a novel "canonical" averaging scheme for the analysis of structure factors of turbulent 
flows. The strategy is to determine the scale r c at which the dissipation rate is the largest in a given 
turn-over time series. This specific scale r c translates into a specific "phase" in the logarithm of the 
scale which, when used as the origin, allows one to phase up the different measurements of a structure 
factor S p (r) = A v (er) p/3 in different turn-over time realizations. We expect, as in Laplacian growth 
and in rupture, that the log-periodic oscillations will be reinforced by this canonical averaging. 
Demonstrating unambiguously the presence of log-periodicity and thus of discrete scale invariance 
(DSI) in turbulent time-series would provide an important step towards a direct demonstration of 
the Kolmogorov cascade or at least of its hierarchical imprint. 



to appear in the Proceedings of the Seventh European Turbulence Conference (ETC-7), June 30- July 3 (1998) (Pub- 
lished by Kluwer, U. Frisch, editor) 
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Scaling laws provide a quantification of the complex- 
ity of turbulent flows. For instance, the second order 
longitudinal structure function has the form ((v r ) 2 ) — 
Cff(er) 2 / 3 , according to Kolmogorov 1941 |jj where r, 
which lies in the inertial range, is the scale at which ve- 
locity differences are measured, and e is the mean rate of 
energy dissipation per unit mass. Dimensional analysis 
shows that ((v r ) 2 ) = (er) 2 / 3 F(!Ji, r/L), where F(x, y) is a 
universal function to be determined, L is the external or 
integral scale. Kolmogorov's assumption is that, for the 
Reynolds number -ft — > oo and r/L — > 0, F(x,y) goes to 
a constant Ck- This is the so-called complete similarity 
of the first kind Q with respect to the variables 3? and 
r/L. 

The existence of the limit of F(8t,r/L — > 0) has first 
been questioned by L.D. Landau and A.M. Obukhov, 
on the basis of the existence of intermittency - fluc- 
tuations of the energy dissipation rate about its mean 
value e. Indeed, Barenblatt's classification leads to the 
possibility of an incomplete similarity in the variable 
r/L. This would require the absence of a finite limit 
for F($t, r/L) as r/L — > 0, and leads in the simplest case 
to the form ((v r ) 2 ) = Cif(er) 2//3 (f)", where a is the 
so-called intermittency exponent, believed to be small 
and positive. If a is real, this corresponds to a simi- 
larity of the second kind 0. Incomplete self-similarity 
HQ may stem from a possible Sft-dependence of the 
exponents. The case where a is complex, leading to 
(( Vr ) 2 ) = Cif(er) 2 / 3 coB[o/log(r/L)], could be 

termed a similarity of the third kind, characterized by 



the absence of limit for F(3t,r/L) and accelerated (log- 
periodic) oscillations. To our knowledge, Novikov has 
been the first to point out in 1966 that structure functions 
in turbulence should contain log-periodic oscillations [m. 
His argument was that if an unstable eddy in a turbulent 
flow typically breaks up into two or three smaller eddies, 
but not into 10 or 20 eddies, then one can suspect the 
existence of a prefered scale factor, hence the log-periodic 
oscillations. They have been repeatedly observed but do 
not seem to be stable and depend on the nature of the 
global geometry of the flow and recirculation |]|| as well 
as the analyzing procedure. 

The theory of complex exponents and log-periodicity 
has advanced significantly jjj in the last few years. Com- 
plex exponents reflect a discrete scale invariance (DSI), 
i.e. the fact that dilational symmetry occurs only un- 
der magnification under special factors, which are arbi- 
trary powers A™ of a prefered scaling ratio A. Complex 
exponents have been studied in the eighties in relation 
to various problems of physics embedded in hierarchi- 
cal systems. In the context of turbulence, shell mod- 
els construct explicitely a discrete scale invariant set of 
equations whose solutions are marred by unwanted log- 
periodicities. Only recently has it been realized that dis- 
crete scale invariance and its associated complex expo- 
nents may appear "spontaneously" in euclidean systems, 
i.e. without the need for a pre-existing hierarchy. Sys- 
tems that have been found to exhibit self-organized DSI 
are Laplacian growth models |q|, rupture in heteroge- 
neous systems , earthquakes [|lO| , animals fll|| (a gen- 
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eralization of percolation) among many other systems. In 
addition, general field theoretical arguments |ll]] indicate 
that complex exponents are to be expected generically for 
out-of-equilibrium and/or quenched disordered systems. 
This together with Novikov's argument suggest to revisit 
log-periodicity in turbulent signals. Demonstrating un- 
ambiguously the presence of log-periodicity and thus of 
DSI in turbulent time-series would provide an important 
step towards a direct demonstration of the Kolmogorov 
cascade or at least of its hierarchical imprint. For partial 
indications of log-periodicity in turbulent data, we refer 
the reader to fig. 5.1 p.58 and fig. 8.6 p.l28ofRef. §, 
fig.3.16 p. 76 of Ref. fig.lb of Ref. || and fig. 2b 
of Ref. §J. 

It is a common observation that the oscillations, if 
any, "move" when changing the length of the signal over 
which the averaging is carried out. They thus have the 
aspect of noise. However, previous numerical simulations 
on Laplacian growth models |8|,|l5|| and renormalization 
group calculations jll) have taught us that the presence 
of noise modifies the phase in the log-periodic oscillations 
in a sample specific way leading to a "destructive inter- 
ference" upon averaging. In the turbulence context, we 
propose that one realization corresponds approximately 
to a signal measured over one turn-over time scale L/vjj. 
In contrast to this sample specific phase dependence, we 
stress that the prefered scaling ratio A has universal prop- 
erties. 

It is thus important to carry out an analysis on each 
sample realization separately, without averaging, as has 
been demonstrated to work for other systems || . An en- 
ticing alternative is to introduce a new averaging scheme 
that does not destroy the oscillations. The standard aver- 
aging procedure, that we could term "Grand canonical" 
Jl6|, is known to introduce spurious sample-to-sample 
fluctuations of relative amplitude proportional to L~ d l 2 
in d dimensions. In contrast, the concept of "canonical" 
averaging [jl6[ consists in identifying, for each realization, 
the corresponding specific value of the critical control 
parameter K^. The natural control parameter then be- 
comes A = (K — K^)/K^ and the averaging can then 
be performed over the different samples with the same 
A. This should then lead in principle to a "rephasing" of 
the log-oscillations. This "canonical averaging" has been 
demonstrated for log-periodic signatures of the acoustic 
emission precursors prior to rupture and in Laplacian 
growth models Eq| . 

We propose to adapt this "canonical" averaging 
scheme to the analysis of structure functions of turbu- 
lent flows. There are probably several possible schemes 
to implement it. Let us suggest here one based on the 
energy dissipation rate. The strategy is to look for a ref- 
erence quantity that is specific to a given turn-over time 
realization. For critical phenomena, a natural candidate 
is the susceptibility whose maximum determines the sam- 
ple specific critical point location ■ For turbu- 



lence, we suggest to determine the scale r c at which the 
dissipation rate is the largest in a given turn-over time 
series. This can be derived by a direct measurement of 
the velocity gradient at small scales or from the third- 
order structure function Ss(r) = — |er, which obeys the 
exact four-fifth Kolmogorov law under a set of assump- 
tions jD. This specific scale r c translates into a specific 
"phase" in the logarithm of the scale which, when used 
as the origin, allows one to phase up the different mea- 
surements of a structure function S p (r) — A p (er) p ^ 3 in 
different turn-over time realizations. We expect, as in 
Laplacian growth models and in rupture, that the log- 
periodic oscillations will be reinforced by this canonical 
averaging. 

What could be the mechanism that creates these char- 
acteristic scales? There are undoubtly the integral length 
L and the scales associated to the dissipation range. But 
what could produce an approximate geometrical series of 
scales in the inertial range? We conjecture two possible 
routes. The first one is inspired from a recent discovery 
that the continuous nonlinear Einstein partial differen- 
tial equations of general relativity in the presence of a 
scalar field self-interacting through gravitation may gen- 
erate a log-periodic spectrum of black hole masses with 
develop according to a log-periodic self-similar time dy- 
namics Jl7[ . The mechanism might result from the ex- 
istence of a limit cycle in the renormalization group de- 
scription of a field close to the negative density limit (in 
turbulence, could this be obtained from a negative effec- 
tive viscosity?). The other route is that scale invariant 
equations that present an instability at finite wavevec- 
tor k decreasing with the field amplitude may generate 
naturally a spectrum of internal scales. An example is 
^ = — 2v|gi — w 2 fjr- This equation is scale invariant 
in the sense that if v(t, x) is a solution, then j 2 v(t, jx) is 
also a solution for arbitrary 7. A linear stability analysis 
shows that a mode Uoe^'e grows with a = 2vok 2 —v 2 ,k 4 , 
i.e. the most unstable mode occurs at finite k m „ = — 3=. 
Thus, a finite characteristic scale appears that is com- 
pletely controlled by the amplitude of the field. Starting 
from an approximate homogeneous level vo, the insta- 
bility produces a large scale fc 2?r = 2n^vo. As dips 
in the field develop, the amplitude there decreases and 
the corresponding instabilities will create smaller length 
scales, and so on. Preliminary simulations Jl8[ confirm 
this intuitive picture : the resulting DSI field is seen to 
result from a cascade of instabilities with characteristic 
wavelengths controlled by the amplitude. 

I hope that these conjectural ideas will stimulate fur- 
ther works on these fascinating log-periodic structures in 
turbulent signals. I am grateful to U. Frisch, L. Gil, N. 
Goldenfcld, A. Johansen, A. Noullcz and G. Simms for 
discussions. 
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